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^ I Abstract 

^ ' A statistical theory of cholesteric liquid crystals composed of short 

pH . rigid biaxial molecules is presented. It is derived in the thermody- 

O ! namic limit at a small density and a small twist. The uniaxial (biax- 

ial) cholesteric phase is regarded as a distorted form of the uniaxial 
(biaxial) nematic phase. The chirality of the interactions and the 
^ . implementation of the inversion to the rotation matrix elements are 

i discussed in detail. General microscopic expressions for the elastic 

■ constants are derived. The expressions involve the one-particle dis- 

. tribution function and the potential energy of two-body short-range 

interactions. It is shown that the elastic constants determine the twist 
of the phase. The stability condition for the cholesteric and nematic 
^ ■ phases is presented. 

The theory is used to study unary and binary systems. The tem- 
perature and concentration dependence of the order parameters, the 
^ ■ elastic constants and the twist of the phase are obtained. The possi- 

bility of phase separation is not investigated. 

Key words: Liquid Crystals; Cholesterics; Elastic Constants; Mix- 
tures. 



1 Introduction 

The cholesteric phase can be considered as a special case of the nematic 
phase [1]. The long axes of the anisotropic molecules are on the average 
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aligned parallel to each other within planes. The direction of this alignment 
rotates smoothly as one proceeds in a direction perpendicular to the parallel 
planes. Such a phase (the twisted uniaxial nematic phase) will be called the 
uniaxial cholesteric phase. By analogy, the twisted biaxial nematic phase 
will be called the biaxial cholesteric phase. We note that in the case of the 
twisted biaxial nematic phase three twists around three perpendicular axes 
of biaxial symmetry should be considered simultaneously. 

There are many models of the uniaxial cholesteric phase where molecules 
are assumed to be uniaxial. In 1970 Goossens [2J showed within the extended 
Maier-Saupe model that the quadrupole interactions give rise to the twist. 
In 1976 Straley [5] presented a statistical-mechanical theory of the elastic 
constants and of the spontaneous twisting of a cholesteric. He expressed the 
pitch by means of the elastic constants. We note that he gave the geometric 
interpretation of the chiral term which is often present in microscopic theories 
of cholesterics. This term can approximately describe the minimum approach 
distance of a pair of threaded rods. We will derive similar terms in the case of 
biaxial molecules. In 1977 Lin-Liu et al. [1] presented a molecular theory of 
cholesteric liquid crystals and showed the general form of the chiral potential 
energy of interactions for uniaxial molecules. They discussed different types 
of the temperature dependence of the pitch. 

Real molecules forming liquid crystalline phases are never uniaxial and it 
is more realistic to assume less symmetric biaxial molecules. Some properties 
of the cholesteric phase formed by biaxial molecules were investigated in the 
past. In 1974 Priest and Lubensky [5J found the order of the biaxial order 
parameter and predicted a fluctuation instability. Van der Meer and Vertogen 
[6] discussed the dependence of the pitch on the biaxial order parameter. 
They also derived expressions for the temperature dependence of the pitch 
and the elastic constants [T] in the case of the uniaxial cholesteric phases. 
Finally, in 1992 Evans [8j presented a hard body model for chiral nematic 
liquid crystals. The density functional theory was used, and molecules were 
represented by a hard convex twisted ellipsoidal core, with and without an 
encircling isotropic square well. The pitch was found to be density and 
temperature independent with values in the visible region of the spectrum. 
Long range potential softness could account for the increase of the pitch with 
decreasing temperature. 

The phase behaviour of liquid crystalline mixtures has been studied, us- 
ing a number of theoretical methods. In 1980 Sivardiere [9] introduced the 
Ising-like model and obtained a large variety of phase diagrams. Brochard et 
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al. [TU] considered the Maier-Saupe model and gave a cataloque of allowed 
diagrams for mixtures of nematogens. As far as the cholesteric mixtures are 
concerned, the theoretical and experimental works describe usually uniax- 
ial cholesteric phases composed of uniaxial molecules. In 1971 Nakagiri et 
al. [11] studied the helical twisting power in nematic-cholesteric mixtures 
as a function of concentration. It was found that, as the concentration of 
the cholesteric material increases, the twisting power of the mixtures in- 
creases lineary up to a certain concentration, beyound which it increases 
more slowly, taking a maximum value at some specific concentration. Then, 
in the high-concentration region, the twisting power falls. We note that typ- 
ically the additivity rule is satisfied as described in Sec. O In 1977 Lin-Liu 
et al. ^ presented a molecular theory of binary cholesteric mixtures. They 
considered uniaxial molecules and derived a formula showing how the pitch 
depends on temperature and concentration. Under certain conditions the 
formula reduces at fixed T — Tc to a simple quadratic rational fraction in the 
composition. 

The forming of the cholesteric phase is closely connected with chirality. 
Some authors tried to establish quantitative relationships between molecu- 
lar properties and measurable properties that result from chiral molecular 
structures. Osipov et al. [12] proposed an intrinsic molecular chirality ten- 
sor based only on nuclear position. The chirality tensor gives rise to two 
universal chirality indices, the first givig information about absolute chiral- 
ity, and the second about anisotropy of the chirality. Recently Harris et al. 
[13\ showed that any chiral measure of a geometric object is a pseudoscalar 
and must involve three-point correlations that only come into play when the 
molecule has at least four atoms. In general, a molecule is characterized by 
an infinite set of chiral parameters. However, one can also consider chirality 
of interactions [IH [T5] and we will use this approach. 

Our aim is to describe the uniaxial and biaxial cholesteric phases that 
consist of biaxial or uniaxial molecules. We would like to derive the micro- 
scopic expressions for the elastic constants and the pitch. We will investigate 
the inversion in the context of chirality. Our paper is organized as follows. In 
Sec. m we present a phenomenological continuum theory of cholesteric liquid 
crystals. In Sec. [3] we describe a statistical theory of cholesteric phases that 
consist of rigid biaxial molecules, and in Sec. IHwe derive general expressions 
for the elastic constants and the cholesteric pitch. A condition of stability is 
obtained that concerns both nematics and cholesterics. Exemplary calcula- 
tions are presented in Sec. 0, were the Corner potential energy is applied. By 
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means of symmetry considerations we will identyfy main cliiral terms for bi- 
axial molecules. In the limit of uniaxial molecules we will recover the results 
from P]. In Sec. [6] we summarize the results of this work. 



2 Phenomenological approach 

In this section we will describe a uniform phase from a phenomenological 
point of view [16]. We assume that at every point r inside a considered 
phase we can define three orthonormal vectors {L{r), M{r), N{r)) reflecting 
orientational properties of this phase. In the case of the biaxial phase they 
determine directions of its two-fold axes of symmetry. The vectors (L, M , N) 
create the local frame which can be expressed by means of a space-fixed 
reference frame (e^, e^, e^) as 

L = L^e^, M = M„e;, N = N^e^, (1) 

where repeated indices imply summation. The completely ordered uniform 
phase is described by (L, M, N) = {e^, Cy, e^). 

Let us call Fa the free energy due to the distortion of the local frame 
(L, M, A^). A general form of its density was derived in [16] in the case 
of small distortions. When a considered phase has a D2 symmetry group 
(the biaxial cholesteric phase) we get 

fd = KuDn + K22D22 + Ks3Ds3 

+^Kini{Dn)^ + \kui2{Di2)^ + ^Kr^uiDrs)^ 

+ \k212i{D2i)^ + \K2222{D22f + ^^^2323 (^^23)' 

+ ^7^3131 (/^3l)' + ii^3232(/^32)' + ^^3333(^^33)' 
1 

+ -(^1221 + ^1122) (^12^21 + ^11^22) 
+ ^(i^l331 + i^ll33) (^13^31 + DiiD-iz) 
+ -(-^2332 + -^2233)(-D23-D32 + -D22-D33) 
+ J(i^l221 - ^1122) (1^12^1 - ^111^22) 
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+ 2 (-^1331 ~ -f^ll33)(-Dl3i?3l 



^11^33) 




-D22-D33) 



(2) 



where Kij, Kijki = Kkuj, Lijk = Ljik are the elastic constants, 



Dii = LaMf^daNp, D12 = LaNpdaLp, Di3 = LaL/sdaM/s, 

Dsi = N^MpdaNp, Ds2 = NaNpdaLp, Dss = NaLpdaMp. (3) 



The terms with Ku give 3 bulk terms, the terms with Lijk give 3 surface 
terms, the terms with Kij^i give 12 bulk and 3 surface terms of the form 



The total numbers of bulk and surface terms are 15 and 6, respectively. 

If a considered phase possesses a symmetry group (the uniaxial 
cholesteric phase), the number of the elastic constants is smaller, because 
some constants from the previous case become dependent or zero. Let the z 
axis be oriented along the axis of symmetry. Then the distortion free-energy 
density has the form of the Frank expression 



da{LpdpLa — LadpLp) 
da{NpdpNa-N^dpNp) 



2{D23D32-D22Dss), 
2(L>3lA3- ^11^33), 

2{Di2D2i-DiiD22). 



(4) 



+ ^K2[N ■ (V X iV)]2 + ^KslN X (V X N)]^ 




(5) 



where the relations among the nonzero elastic constants are 
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1 — 



1212 — 



2121; 



K2 — -ft'llll — K22221 

K3 = Ksisi = K3232, 




(6) 



Therefore in the case of the uniaxial phase we have 4 bulk (from Kq to K^) 
and 2 surface terms {K4 and K^). 

In the continuum theory of uniaxial nematic liquid crystals three basic 
types of deformations i.e. splay, twist and bend, appear, which extract from 
the distortion free energy terms with Ki, K2 and K^, respectively. Thus each 
constant Ki must be positive; if not, the undistorted nematic conformation 
would not correspond to a minimum of the free energy Fd- In [T7] 18 basic 
deformations proper for the continuum theory of biaxial nematics were given. 
They were divided into five groups and connected with relevant elastic con- 
stants: 3 twists (for Kiiii), 6 splays and bends (for Kijij), 3 modified twists 
(for Lijk) and two groups of 3 double twists (for Kajj and for Kyji). Inside 
the formulas for deformations a parameter q was used (1/g was a certain 
length). Small q meant a small deformation and a conformation close to the 
uniform one (L(°), M(°), iV(°)). The vectors of the local frame were expanded 
into a power series with respect to q: 



It appeared that the most important in Eq. ([7]) were terms linear in q. They 
were sufficient to calculate the distortion free energy up to the second order 
in q, and to calculate the elastic constants of biaxial nematic liquid crystals. 
We add that for Ku one should use the same deformations as for Kuu because 
the relevant terms have different dependence on the parameter q and they 
do not mix. 

3 Microscopic approach 

In this section we focus on the microscopic analysis of cholesteric liquid crys- 
tals. Let us consider a dilute gas of N molecules contained in a volume V, at 



M 



N 



L 



1(0) ^ ^IH) ^ ^2£(2) ^ _ _ _ ^ 

M (0) + qM^i) + ^2^(2) ^ _ _ _ ^ 



(7) 
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the temperature T (in Kelvine) . Let it be a binary mixture of molecules 
A and Nb molecules B, where = A^^ + Nb- We assume that molecules are 
rigid blocks with three translational and three rotational degrees of freedom. 
The state of a molecule i is described by a vector of the position rj and the 
orientation Ri = {(f)i,6i,tpi), where (pi, 6i and ipi are the three Euler angles. 
On the other hand, one can use the set of the three orthonormal vectors 
(/, m, n). In a space-fixed reference frame (e^., e^, e^) we can express them as 

Let m/, J/a;, Jiy and Jiz denote the mass of a molecule I {I = A oi B) 
and the three moments of inertia, respectively. Apart from that, we denote 
■u = r 2 ~ = wA, Ru = {(f), 9, ip) are the three Euler angles for A. In fact, 
ip is not used, because A = (sin6'cos(/), sin6'sin(/),cos^) and Ru will appear 
inside the rotation matrix elements -D^o- 

We assume that molecules interact via two-body short-range forces which 
depend on the distance between molecules and their orientations; ^"(2 (m, Ru, Ri, R2) 
gives the potential energy of interactions (/, J = A, B). Now we would like 
to discuss in detail symmetries of the energy ${2 molecules, because it 
is crucial for the forming of the cholesteric phase. We note that the Euler 
angles always enter any formulas via the standard rotation matrix elements 
Djll There is a problem how to incorporate the inversion Ci into this for- 
malism because the inversion can not be expressed by means of rotations. 
However, one can try to simulate the inversion. We start from the vector A. 
After the inversion operation it gives —A. Let us denote the Euler angles 
after the inversion by IRu- 

= D^li(t> + n,n-e, = (9) 

We used an unknown function \l/ of the angle ip because the third Euler angle 
is not determined before and after the inversion. Note that the relation 
is often generalized to the form [IB] 

DililRu) = {~l)'D%{Ru)- (10) 

The relation fllUp in general is not true because it cancels the matrix elements 
D^^j^ with / odd. On the other hand such elements have to be used as described 
in dH. 
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Let us rewrite the relation in the form 



(11) 



Thus, for the matrix elements D^j^q the inversion can be connected with the 
rotation Rz{tt) = (0, tt, 0) and the C2 symmetry group. This means that there 
is no difference between D^o and -Dooh within the considered formalism. 

We postulate that the invariance with respect to the inversion means the 
invariance with respect to the three separate rotations Rxi'n'), Ryi'^) and 
R,{tt) and one should replace D^^^lilRu) with Dll^RuRxij^)) , Dll^RuRyiji)) 
and D^^I{RuRz{tt)), respectively. This connects Ci with the D2 symmetry 
group . This means that there is no difference between D2 and D2h within 
the considered formalism. It is known that Rxi"^) = Ry{T^)Rz{'^) so only two 
rotations have to be used. 

Now we are in a position to discuss all symmetries of the potential energy 

1. Translational invariance: It is satisfied because $(2 depends on u. 

2. Rotational invariance: ${2 should not depend on a choice of a reference 
frame. It means that for any rotation R 



We add that generally it does not have to be satisfied because the 
interactions between two molecules can be modified by the presence of 
others molecules, especially in an ordered phase. 

3. Invariance with respect to the permutation of identical molecules: 



4. Invariance with respect to the symmetry operations of molecules: for 
biaxial molecules we apply operations from the symmetry group. 
For the molecule / 



^[i{u, RRu, RRi, RR2) = ni{u, Ru, R2). 



(12) 



IRu, R2, Ri) = Kiu, Ru, Ri, R2). 



(13) 




{u,Ru,RiRz{n),R2) 




{u, Ru, Ri, R2) 



(14) 




{u,Ru,RiRy{n),R2) 




{u, Ru, Ri, R2), 



(15) 
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and for the molecule J 

^{i{u, Ru, Ri, R2Rzi7r)) = $ R^, Ri, R2), (16) 

^{'^{U, Ru, Rl, R2Ry{TT)) = (m, Ru, Ri, R2). (17) 

If the molecule I is uniaxial, we should add operations from the Dooh 
symmetry group. For any a 

^{i{u, Ru, RMa), R2) = ^iUu, Ru, Rl, R2). (18) 

5. Invariance with respect to complex conjugation: $(2 should be a real 
function. 

6. Chirality of the interactions: the interactions between molecules are 
nonchiral when 

¥^i{u, IRu, IRi, IR2) = $(2 (m, Ru, Rl, R2), (19) 

where \R denotes the Euler angles R after the inversion. In the opposite 
case the interactions are chiral and this leads to the forming of the 
cholesteric phase. 

The microscopic free energy of the binary mixture has the form [19] 
j3F = /f/(l)G'/(l){ln[G,(l)A/]-l} 

I=A,B-' 

-\ E ld{l)d{2)G:{l)Gj{2)fii, (20) 

^ I,J=A,B-' 

where Gi{l) = Giifi, Ri) are the one-particle distribution functions with the 
normalizations 

j d{l)Gj{l) = Nj, (21) 

d{l) = dridRi, = exp(— /S^jg) ~ 1 the Mayer functions, P = l/ksT, 
and ^ 

A7=(^^) {m]Ji,JiyJi,y''\ (22) 

Our set of state variables consists of T, V , Na, and A''^. The free energy f l20l) 
consists of the ideal terms (with A/) and the excess terms directly related to 
intermolecular forces. The ideal terms are that of the ideal gas. 
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The expression fl20p was derived systematically for binary mixtures from 
the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy equations in the ther- 
modynamic limit {N — »• oo,V — *• oo,N/V = const) |T9]. The two-particle 
distribution functions were expressed in terms of the one-particle distribu- 
tion functions and the two-particle correlation functions of the simple form 
exp(— /3$(2)- This assumption guarantees the proper limit of the unary sys- 
tem. 

The equilibrium distributions Gj minimizing the free energy fl20l) satisfy 
ln[G',(l)A/] - E / f^(2)Gj(2)AY = const. (23) 

J=A,B-' 

In the homogeneous phase the distribution function Gi do not depend on 
the position of a molecule and Gi{l) = Goi{Ri). In order to obtain Gqi one 
should solve the equations fl23|) together with fl2T|) . 

In order to define the microscopic distortion free-energy density fd one 
should also identify the homogeneous free-energy density /q. Note that the 
terms with A/ in fl2Ul) are local: the integrands involve the distribution func- 
tion for a single point only. Other terms couples the distributions at neigh- 
bouring points. The dependence of the free energy on the spatial variations 
of the ordering will be found by expanding this terms in the gradients of Gj 
[3]. Substituting the Taylor expansion of Gj{ri + u) into (120|) we get 

F = j cir[/o(f) + /,(r)], (24) 

where 

/5/o(r1 = E I dRGi{r,R){\n[Gi{r,R)Aj]-l} 

I=A,B-' 

-\ Yl fdRidR2duGj{f,R^)Gj{f,R2)f(l (25) 

^ I,J=A,B-' 

f3fd{^ = E j dR^dR2duGi{r,Ri)[{u-V)Gj{r,R2)]fii 

+ 7 E / dRidR2du[{u ■ V)Gi{r, R,)] 

^ I,J=A,B'' 

x[{u-W)Gj{r,R2)]fii. (26) 

Note that an integration by parts has been used to combine the second or- 
der terms and the surface terms have been neglected thanks to the ther- 
modynamic limit. The definition ( l25l) is equivalent to that by Poniewierski 
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and Stecki [20]. This is a well-founded assumption if we also assume slow 
variations of the vectors {L,M,N). We will also restrict the one-particle 
distribution function Gj to the class of Gqj functions. This method was suc- 
cesfuUy used in the past [211 [22] • As we expect, for the homogeneous phase 
fd becomes equal to zero. 

It was shown in [231 [17] that in the case of the homogeneous biaxial 
nematic phase composed of biaxial molecules, the one-particle distribution 
function Gqi depends on four arguments: 

Goi{R) = Goi{l -e^J- e;, n-e^,n- e^). (27) 

We postulate that the distribution of the distorted phase Gi{f,R) can be 
written as 

Gi{f, R) = Goi{l- L(f), f - N{f),n ■ L{r),n ■ iV(f)), (28) 
where the reference frame (e^, Cy, Cz) is replaced with the local frame {L{r), M[r), N{f)). 

4 Elastic constants 

Now we are in a position to substitute the basic deformations into the mi- 
croscopic distortion free-energy density ([2^ and to the phenomenological 
distortion free-energy density ([2|) and (jSj). As a result of the comparison we 
get the microscopic formulas for the elastic constants. To make them more 
compact we write 

= diGoila + dsGoin^, 
Wi = d2Goila + d,Goina. (29) 

Microscopic expressions for the chiral elastic constants of the biaxial choles- 
teric phase follows: 

f3Kn = -\l dR.dR^duu, ^ f(,'Goi{Ri)W,'y, 

^ •' I,J=A,B 

PK22 = -\ 1 dRrdR2duUy fliGoi{Ri){Ui-Wi), 

^ '' I,J=A,B 

(3Kss = \ I dR,dR2duUz J2 fuGoi{Ri)U^y. (30) 

I,J=A,B 
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Microscopic expressions for the nonchiral elastic constants are the following. 
The first group is, 

rdR^duul fuK^^ (31) 

I,J=A,B 



PKnn = 


M 


j3K2222 = 






II 



I,J=A,B 



cIJttI ttJ 



(33) 

I,J=A,B 



The second group is, 

/?i^i2i2 = UdR.dR^duul Y: fL'iK-WDiK-Wi), (34) 

^ '' I,J=A,B 

/3i^i3i3 = ^ / dR,dR2duul Y fHKK^ (35) 

^ ^ I,J=A,B 

/3i^2i2i = ^ / dR,dR2duul Yl fuWlyW^y, (36) 

^ I,J=A,B 

PK2S23 ^ \ j dR^dR^duul Y fiiUiyK^ (37) 

^ I,J=A,B 

/3i^3i3i = ^ / dR^dR2duul Y fii^ly^i^ (38) 
/9i^3232 = IfdRMduul Y fl2{Ul-WL){Ui-Wi). (39) 

^ I,J=A,B 

The third group is, 

-^123 = -^231 — -^312 = 0. (40) 

The fourth group is, 

(5Kii22 = dRidR2duUxUy 

X E fl2mL-Wi,)Wi + Wl^{Ui-W2i)], (41) 

PK2233 = ^ J dRidR2duUyUz 

X E fii[-uiy{ui-wi)-{ui-wl)ui^i (42) 

12 



PKu33 = - dR,dR,duu^u, Yl fiii-UiyWi-^lyKy (43) 

^ I,J=A,B 

The fifth group is, 

K122I = Kii22i i^l331 = -^^1133) -f^2332 = -^^2233- (44) 

In the case of the uniaxial cholesteric phase the expression for the chiral 
elastic constant is 

mo = \l dRrdR^duUy fiiGoi{Ri)Wi. (45) 

^ ■' I,J=A,B 

The expressions for the nonchiral elastic constants are as follows: 





- J dRidR2duul 


I,J=A,B 


(46) 


PK2 = 


- J dRidR2duUy 


E f(2wLw2i, 

I,J=A,B 


(47) 




- J dRidR2duul 


E fl2wLw2i, 

I,J=A,B 


(48) 




l(Ki + K2), 




(49) 


= 


0. 




(50) 



Let us show how the elastic constants determine the state of cholesteric liquid 
crystals. It is very important that, when we look for the physical state of 
the liquid crystal in the thermodynamic limit, we should remove all surface 
terms from the phenomenological free-energy densities. The reason is that 
in the thermodynamic limit all surface terms are negligible. Wc used the 
surface terms to derive the expressions for the elastic constants to assure the 
consistency because the surface terms are present in the hidden form in the 
microscopic free-energy density. The equilibrium distortion of the uniaxial 
cholesteric phase is a pure twist 

N{f) = [0, — sm{qx),cos{qx)] (51) 

where q = Kq/K2, P = 27r/|g| is the cholesteric pitch and = —Kq/2K2- 
Note that the sign of Kq distinguishes between right- and left-handed helices. 
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For the biaxial cholesteric liquid crystals the equilibrium distortion is 
composed of three twists with respect to the orthogonal axes. In the limit of 
small distortions the vectors of the local frame have the form 



L = {l,q3Z,-q2y), 

— * 

M = {-q3Z,l,qix), 

— * 

N = {q2y,-qix,l), 



where qi 



Wi/W, 

w 



— -f^llll-f^2222-f^3333 + 2i^li22^1133-^2233 

— -f^llll-f^2233 ~ -f^2222-f^ii33 — -f^3333-f^ll22) 
Wi — i^ll(i^^2222-f^3333 ~ -^2233) 

+ -f^22(-f^l 133-^^2233 ~ -f^ll22-f^3333) 
+ -f^33(-^1122-f^2233 ~ -f^l 133 -^2222 ) , 
W2 — -f^'ll(-f^ll33-f^2233 ~ -^^1122-^3333) 
+ -ft'22(-ft'llll-ft'3333 - -^1133) 
+ -^33 (-^1122-^1133 ~ -^1111-^2233)) 
Ws — -ft'll(-ftril22-ft^2233 — -^1133-^2222) 

+ -f^22(-f^ll33-f^ll22 — -f^l 11 1-^^2233) 



+ -f^33(-f^l 111 -^2222 ~ -^1122)- 

The distortion free-energy density at the minimum is 

fd — [2-ft'ii-ft^22 (-^3333-^1122 — -^1133-^2233) 

+ 2A'iiA'33(A'2222-^1133 ~ -^1122-^2233) 
+ 2^22-^^33 (-^1111-^2233 " -^1122-^^1133) 
-^2222-^3333) 
-^llll-R'3333) 



+ K^^ (-^2233 



133 



(52) 



(53) 



+ KUKI,22 - A'lllli^2222)]/2Py^ 

It is important that a global minimum of fd exists only if the matrix 



(54) 



i^iiii 

-f^ll22 
-f^ll33 



-f^ll22 
-f^2222 
-f^2233 



-f^ll33 
-f^2233 
-f^3333 



(55) 



14 



is positive definite. In the opposite case there is no stable biaxial cholesteric 
(and nematic) phase. Taking the uniaxial phase limit, we get simpler con- 
ditions. We can say that there is a stable cholesteric or nematic phase only 
if 

Ki < 3K2. (56) 

As far as we know almost all theories and all real and computer experiments 
are in the agreement with the relation (1561) . The equality Ki = 3K2 appears 
in some theories with hard molecules f21\ . 



5 Exemplary calculations 

The aim of this section is to express the elastic constants by means of the 
order parameters which can be measured in experiments. We will apply 
the Corner potential energy of the interactions because in principle it al- 
lows detailed calculations without any additional approximations. On the 
other hand, it is quite realistic. The Corner potential energy has the form 
^{2{u/a^'^), where a^"^ depends on orientations -Ri, R2 and For a^'^ one 
can write the general expansion proposed by Blum and Torruela ^18j. It in- 
volvs the 3-j Wigner symbols and the standard rotation matrix elements. The 
same expression was used to describe the interactions of biaxial molecules in 
|23j . In the case of biaxial molecules the lowest order terms of the expansion 
give 

+ ai({A ■ hf + aii{A ■ l2f + cji\h ■ k? 

+ aliih ■ n2f + (Tli{ni ■ fs)^ + (jl\ni ■ n2){ni x ^2) ■ A 

+ (Tj([{li ■ n2)(/i X n2) ■ A - (mi ■ n2)(mi x 712) ■ A] 

+ ^liii.^! ■ k){ni XI2) ■ A - {ni ■ m2)(^i x m2) ■ A] 

+ cr8-^[(/i ■ hWi X I2) ■ A + (mi ■ fh2){rrii x m2) ■ A 

- (fi ■ m2)(ri X m2) • A - (mi ■ r2)(mi x I2) ■ A]. (57) 

There are 13 molecular parameters that determine the main features of the 
interactions between biaxial molecules. We would like to add some comments 
to this long expression. The terms with ctq"', af/, (7(2 5 ^"^^ describe the 
nonchiral interactions of uniaxial molecules. The term with describes 
chiral part of their interactions. A similar term was used in the past [U ■ 
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The terms with cTq'^, erf/, erf/, crg'^, cr|/, cr|/, cr|'^, cr|/, and cXg/ describe the 
nonchiral interactions of biaxial molecules [231 [IZj. The additional terms 
with (Tg"^, erf/, erf/, and o"|'^ allow to decribe the chiral interactions of biaxial 
molecules. As far as we know they were not presented in the literature. 
Note that the terms with A couple the translational and rotational degrees 
of freedom. 

In order to set the proper values of the parameters crf'^ one should in- 
vestigated the equipotential surfaces of the potential energy. The excluded 
volume method described in [23] can be helpful. Molecules I are described 
by 



.11 — ^11 — ^11 ^11 ^11 — ^11 — ^11 ^11 



crf^ ^ aii = ail ai\ af^ = aff = aff, ai' . (58) 

The interactions between molecles / and J we describe by 

= + < = ^V, erf/ = cr/^, 

ai^ = {ai^ + a/^)/2, a^/ = , a|/ = a/^, 
^l' = (c^V + ^/")/2, aM = aM = (af + cri-^)/2, 
^i' = + ^i')/2, erf/ = af/ = (af^ + a/^)/2, 

ai^ = {ai^ + ai-^)/2. (59) 

The parameters cxf are the molecular constants that in our model determine 
the ordering of the phase, the value of the elastic constants and the twist of 
the phase. As far as the functional dependence of the potential energy on 
the u/a^^ is concerned, we have many possibilities and we will mention two 
of them, together with a function Bs{T) defined as 

Bs{T)= dxx'fii{x)= rfa;x'[exp(-/5$f/(x)) - 1]. (60) 
JO Jo 

Let e denotes a depth of the potential energy (we assume for simplicity that 
it is the same for both types of molecules). 

1. The soft-core potential energy 

^{i{ula^-^) = e{cr^'/ur, (61) 

-1 /m-s-l\/ e 
BsiT) = -—T^ — - . (62) 



s + 1 V m J KksT 
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2. The square-well potential energy 

{+00 for iu/a^-^) < 1 
-e for 1 < {u/a^-^) < Rsw (63) 
for {u/a^J) > Rsw, 

^^(^) = JTl [("^P('/^^^) - ^) (^^^ - 1) - 1] ■ (64) 

Let us define dimensionless functions fi{R) = Goj{R)V/Ni with the normal- 
ization 

' dRfj{R) = 1. (65) 



In order to describe systems with biaxial symmetry it is convenient to in- 
troduce the complete set of basic functions (invariants) proper for that kind 
of symmetry. They are closely connected with the rotation matrix elements 
and have the form 

C, \ 2+(5po+<5oi/ 
T.J-^y'^''''^^'Dl^i.-u{R), (66) 

where j is a nonnegative integer, /i, u are even. If j is even, then < fi < j 
and < u < j . If j is odd, then 2 < fi < j and 2 < v < j. The invariants 
are real and orthogonal functions 

/ dRF%\R)F',t\R) = S,,S,,S,.8nyi2j + 1). (67) 



The properties of the invariants were described in [17J. It was shown that all 
FjfJ can be expressed by means of {I ■ e^y, (/ ■ e^)^, (n ■ e^Y, and (n ■ e^)^. 
Let us show the most important invariants with j = 2: 

fS'> = [-l + 3(n-e,)2]/2, 

= [-l + (n-e,)2 + 2(n-e;)']V3/2, 
fS! = [-3 + {n-e,Y + 2{l-e,Y + 2{n-e,Y + A{f-e,f]/2. (68) 

For the first time all four invariants (with different normalization) appeared 
in the paper by Straley [26]. Let us define the order parameters for our 
system as 

{f!^)j^ JdRfj{R)F!^J{R). (69) 
17 



Note that for the completely ordered biaxial phase we have {FjfJ)]: = 5^^^ 

f2) ('2) 

thus the most important are {Fqq )i and (-F22 )/ V in the Straley 

notation). 

Let us rewrite (123!) for the one-particle distribution function of the uniform 
phase. It is identical to one in the theory with nonchiral interactions. We 
define 

K''{RuR2) = j dA{a''/al')\ (70) 

X'' = B^{T){al,'fNj/V. (71) 

The kernel K^^ and the distribution functions // can be expressed in terms 
of the invariants 

K^'{R.Q) = Y.llK\iy'Fl^XR). (72) 
ln//(i?) = EE4^'^^.^(^)- (73) 

j flU 

Note that for a^'^ given by (!57|) both sums (!72|) and (173|) are finite series. 
Thus the distribution function // is fully described by 35 coefficients Sjj'J^ . 
We calculate the coefficients S'^J^ from the equations 

S;£'= E A^^E^i^"(i^^;)^forj>0, (74) 

J=A,B P 

(Fqq)/ = 1 (the normalization condition). (75) 

Now we move to the elastic constants. We can expand the distribution 
function fj in an infinite series with respect to invariants 

MR) = EE(^/i^)/^iiH^)(2j + l)/{8rr'). (76) 

j tiu 

All elastic constants can be written as finite sums of the form 

K.. = T.^''{F^)i{Fl^)jBy!i^, (77) 

K.,M = E v"{F;£)i{F^)jA^lii^, (78) 

where 

= kBTB^{T){ai'YNiNj/V\ (79) 

77" = kBTB^{T){ai-'fNiNj/V\ (80) 

The coefficients K\^y^ , Afk^J^a^ and B\l^jll[^ are polynomials in (af'^/a^'^) 
that can be calculated analytically. 
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5.1 Unary phase of biaxial molecules 

For simplicity in our calculations we took into account only the dependence 
on the order parameters (Fqq^) and {f!^). We assumed the square- well 
potential energy of interactions with Rsw = 2. Note that o"o determines the 
length scale, whereas e determines the energy scale. The elastic constants 
will be expressed in e/cxo and c/uq, temperatures will be given in e/ks- We 
assumed that molecules are similar to ellipsoids with three different axes 
((To/2) X (cTo) X (2(To). We used V/N = Sa^ and our set of the molecular 
parameters is 

(Ji = ao/2, 0-2 = 0-0/2, (T3 = -0-0/4, 

0-4 = 0-5 = 0, 0-6 = 0-7 = 0-8 = 0-0/10. (81) 

In our system, on decreasing the temperature we meet the first-order tran- 
sition to the uniaxial cholesteric phase at T = 1.94 and the second-order 
transition to the biaxial cholesteric phase at T = 0.41. The temperature 
dependence of the order parameters (-foo'') ^^nd (^22^) is presented in Fig. [TJ 
The order parameters were used to calculate the temperature dependence of 
the elastic constants. It appeared that the factors ^ and r] introduce very 
strong variations with temperature and it was more comfortable to present 
the elastic constant divided by those factors as shown in Figs. [2] and [31 The 
corresponding twist parameters q and are shown in Fig. HI 

In the uniaxial cholesteric phase there is one chiral constant (Kq) and 
three nonchiral ones [Ki, K2, and K^). The twist parameter q describing the 
pure twist is decreasing with temperature. For some potential energies one 
can obtain also the increasing temperature dependence, more often present 
in cholesterics. We have also the accidental equality Ki = caused by the 
omission of the order parameters {F^J) with j greater than 2. Typically, for 
rodlike molecules we get the inequalities K2 < Ki < K^. 

In the biaxial cholesteric phase we have three chiral constants (-ft'n, -ft'22, 
and -ft'as) and 12 nonchiral ones. On decreasing the temperature, Kq splits 
into Kii > K22 and K33 becomes nonzero. Ki splits into -ft'1212 > -^2121, 
K3 splits into -ft'3232 > -^3131, K2 splits into -ft'2222 > -f^iiii, -^^1122 becomes a 
new independent constant. Other new independent constant are four positive 
K1133, -fCisis, -R'2323, -^3333 and one negative -ft'2233- They are about two orders 
of magnitude smaller then the previous ones. The twist parameters qt are 
quickly increasing with the decreasing temperature and satisfy inequalities 
qi< q2< qs- 
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5.2 Binary phase of uniaxial molecules 

We mixed the prolate molecules interacting via the soft core potential with 
m = 12. The set of molecular parameters is 

a^^ = ao, = 0.2ao, a^"" = -0.2ao, a^^ = -O.lao, 

(To^^ = ao, af^ = 0.3^0, af^ = -0.2ao, af^ = 0, (82) 

where ctq denotes some molecular length and V/N = 4.471(Tq. The elastic 
constants will be expressed in 10~^e/(jQ and 10~^e/(To, the temperature in 

For simplicity reasons we assume that only one phase is present. It corre- 
sponds to the stable solution with the minimum free energy. The possibility 
of phase separation will not be investigated. We checked the stability against 
perturbations of nematic symmetry and the stability against deformations of 
the phase. 

The molecules A and B in the unary phase form the uniaxial cholesteric 
and nematic phase, respectively. The transition temperatures of unary sys- 
tems from the isotropic to the uniaxial cholesteric or nematic phase are 
Ta = 124 and Tb = 162. The temperature dependence of the order param- 
eters (-P2)/ and the elastic constants for unary systems are shown in Figs. [5] 
and [6l The pitch of the system of molecules A is a falling concave function 
of the temperature. 

In the considered range of temperatures the isotropic or uniaxial cholesteric 
phase was present in the mixture. The transition temperature Tc of the 
binary system from the isotropic to the uniaxial cholesteric phase can be 
described approximately as 

Tc= ^iTi. (83) 

I=A,B 

We performed a detailed analysis of the mixture at the temperature T = 100 
where the uniaxial cholesteric phase was present for all concentrations. The 
order parameters and the elastic constants are continous and concave func- 
tions of the concentration xa- The dependence of the elastic constants on the 
concentration is plotted in Fig. [3 The twist of the mixture is approximately 
described by the equation 

q= ^Kli^ (84) 

I=A,B 
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where qj are related to the unary systems at the same temperature. The 
additivity rule fl84p is satisfied for most mixtures between cholesterol deriva- 
tives yj but there are some exceptions [27]. Note that in general qi may be 
positive or negative. 

6 Conclusions 

In this paper we developed the statistical theory for the uniaxial and biax- 
ial cholesteric phases. Rigid molecules interacting via two-body short-range 
forces were assumed. We derived the microscopic formulas for the elastic 
constants and the pitch of the phase. In order to calculate the values of 
the elastic constants one needs the one-particle distribution function and the 
potential energy of the molecular interactions. It is showed that the elastic 
constants determine the twist of the phase. The obtained stability condition 
is to our knowledge not to be found in the literature. In the case of the 
uniaxial cholesteric or nematic phase it has the simple form Ki < 3K2. It 
seems that it is satisfied for all known substances. 

Our theory was applied to the unary and binary systems of molecules 
similar to ellipsoids with three different axes. Thanks to the Corner po- 
tential energy the elastic constants were expressed as a finite series of the 
order parameters. Apart from this the role of the temperature is more trans- 
parent: the function B2{T) determines the order parameters; the functions 
B^{T) and Bi{T) are closely related to the chiral and nonchiral elastic con- 
stants, respectively. The temperature dependence of the order parameters, 
the elastic constants and the twists in uniaxial and biaxial cholesteric phase 
was obtained. 

The results concernig the nonchiral elastic constants Kij^i are the same as 
in the case of nematic phases described in [17] . We focus on the new results 
on the chiral elastic constants Ku. During the transition from a uniaxial to 
a biaxial cholesteric phase Kq splits into Ku and K22-, whereas -ft^ss becomes 
nonzero. Below a certain temperature they satisfy inequalities Ku < K22 < 
K33. The corresponding twist parameters satisfy similar relations qi < q2 < 
q^. Note that those relations depend on the type of molecular interactions 
and do not have to be always satisfied. 

In the binary mixture it appeared that the transition temperatures and 
the twist of the mixture approximately satisfy the simple additivity rules. 

Our theory is a starting point for further researches. It is desirable to 
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go beyound a low density limit where the Mayer function is replaced with a 
better approximation of the direct correlation function C2. The theory can be 
generalized to the case of more component liquid crystalline mixtures. We are 
also waiting for the experimental data on the biaxial nematic and cholesteric 
phases and the elastic constants as, according to Kini and Chandrasekhar, 
experiments are feasible [2H]- The experiments can verify the theoretical 
predictions and provide new problems. 
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Figure 2: Temperature dependence of the chiral elastic constants Ku/^. In 
the uniaxial cholesteric phase (between T = 1.94 and T — 0.41) we have 
Ko = Kii = K22 and ^^33 = 0. In the biaxial cholesteric phase (below 
T = 0.41) Kq splits into two constants and i^33 becomes nonzero. 
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Figure 3: Temperature dependence of the nonchiral elastic constants Kij^i/r). 
In the uniaxial cholesteric phase (between T = 1.94 and T = 0.41) we have 
Ki = Ki2i2 = -ft^2i2i equal to K3 = K3131 = K3232 and K2 = Knn = 
K2222- We have also i^^ii22 = {K2 — Ki)/2. In the biaxial cholesteric phase 
(below T — 0.41) Ki, K2 and K2, split into two constants, K1122 becomes an 
independent constant, five constants become nonzero (i^uss, -^^2233, -^^1313, 

^^2323 and -fr3333)- 
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Figure 4: Temperature dependence of the twist parameters q and qi. In 
the uniaxial cholesteric phase (between T = 1.94 and T = 0.41) we have 
g > 0. In the biaxial cholesteric phase (below T — 0.41) all are positive 
and qi <q2 < qs- 
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Figure 5: Temperature dependence of the order parameters Sa 
(Cholesteric) and Sb — {P2)b (Nematic) for the unary systems. 
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Figure 6: Temperature dependence of the elastic constants for the unary 
systems. 
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